Voltage-Controlled Surface Magnetization of Itinerant Ferromagnet Nii.xCux 
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We argue that surface magnetization of a metallic ferromagnet can be turned on and off isother- 
mally by an applied voltage. For this, the material's electron subsystem must be close enough to the 
boundary between para- and ferromagnetic regions on the electron density scale. For the 3d series, 
the boundary is between Ni and Cu, which makes their alloy a primary candidate. Using Ginzburg- 
Landau functional, which we build from Nii_xCux empirical properties, ab-initio parameters of Ni 
and Cu, and orbital-free LSDA, we show that the proposed effect is experimentally observable. 

PACS numbers: 85.75.-d, 85.70.Ay 



Further evolution of magneto-electronics [H depends 
highly on the availability of materials, in which local mag- 
netization can be turned on and off isothermally by an 
electric voltage. The hopes to achieve this are mainly 
laid on the dilute magnetic semiconductors (DMS) ^] , in 
which the effect was recently demonstrated experimen- 
tally The voltage-controlled ferromagnetic ordering 
in DMS relies on the virtue of doped semiconductors to 
allow external variation of the free-carrier spatial density 
within the semiconductor depletion layer, typically mea- 
sured in dozens of nanometers. The voltage variation of 
the high electron density in metals is possible only within 
the atomic size Thomas-Fermi (TF) surface layer. As a 
result, the voltage controlled ferromagnetism in a metal 
has not been considered a possibility lately. 

We argue that by capacitively charging a metallic fer- 
romagnet one can drive the surface electron subsystem in 
and out of its ferromagnetic state. At this, the electron 
system of the metal has to be paramagnetic at the de- 
vice operation temperature (room temperature, Tq), but 
close enough to the ferromagnetic state on the temper- 
ature and/or electron density scales. The proximity of 
the ferromagnetic transition will play a twofold role: (i) 
the capacitive change in the electron density is relatively 
small, so that the transition has to be sufficiently close in 
order to reach the ferromagnetic region with reasonable 
voltages, and (ii) due to the critical collective spin cor- 
relations, the spin correlation length grows infinitely as 
one approaches the transition point. Consequently, even 
though the injected carriers are spatially limited to the 
TF-layer, the system must develop a much wider surface 
magnetization profile. In this paper we investigate the 
proposed possibility in Nii_xCux. 

The nature of ferromagnetism lies in the competition 
between the kinetic energy and the exchange interaction. 
The kinetic energy of the spatial quantization tends to 
equalize the numbers of spin up and down electrons by 
shifting the fermionic antisymmetry into the spin sector 
of the many-body wavefunction. In turn, the exchange 
interaction does the opposite, struggling to unbalance the 
up and down spins. For itinerant ferromagnets, - the sub- 
class of materials elemental ferromagnets belong to, - the 
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FIG. 1: (a) Schematics showing the paramagnetic chemical 
potentials of the elemental ferromagnets on the kinetic energy, 
e, and the spatial electron density, po (outer shell electrons 
per atom), scales. The inset is the magnified dashed area, 
which shows the graphic solution for Eqs.(|3]) determining the 
spin density, sq, the paramagnetic and ferromagnetic chemical 
potentials, /ip and /.i, and the spin up and down kinetic Fermi 
energies, e]^^, as functions of po for pure Ni, po = pNi. (b)-(d) 
The experimentally observed properties of Nii.xCuxi (b) zero- 
temperature magnetization vs. po; (c) Curie temperature vs. 
Po; (d) the magnetization vs. temperature, T . 



outcome of this competition can be predicted from the 
Stoner criterion According to the Stoner criterion, 
ferromagnets possess a high density of states (DOS) at 
the (paramagnetic) chemical potential. In the 3d series, 
the high DOS is provided by the 3d band on the back- 
ground of the low and wide 4s and 4p bands (see FigH^). 
Thus, the elements with the chemical potential within 
the 3d band, Fe, Co, and Ni, are ferromagnetic, whereas 
the very next element, Cu, is not. Cu has an extra elec- 
tron per atom beyond Ni, so that on the electron density 
scale the ferro-to-para boundary lies between Ni and Cu. 

Despite their different magnetic properties, Ni and Cu 
are very similar from the band structure point of view. 
In a crystalline state both form the face-centered cubic 
lattice with almost the same lattice constant, a (6.69 vs. 
6.83 a.u.). In result, on alloying, Ni and Cu form sub- 
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stitutional solid solutions at all compositions. Conse- 
quently, Nii_xCux can be though of as an all-the-same 
structure solid, with the equilibrium electron density 
varying linearly with x, pq = p^i + x/a'^. Within this 
picture, x and po are interchangeable. 

Other experimentally observed properties of the al- 
loy 0] can be well approximated as: (see FiglTjD-d) the 
linear dependence of the zero temperature spin density, 
so(pO)O), and the Curie temperature, Tc{pq), on po, and 
the Landau dependence of the spin density, sq{pq,T), on 
the temperature, T: 

2so{po,0)/j = {p3d - PoWipsd - Pq), (la) 
Tc{po)/k = {P3d ~ Po)0{p3d ~ Po): (lb) 
so{po,T)/so{po,Q) = (l-f)i/20(p3d-po), (Ic) 

where T = T/Tc{po), 9 is the Heaviside step- function, 
P3d = PNi -I- 0.53/a'^ is the position of the 3d band edge 
on the density scale, and 7 and k are the slopes of the 
2so(po,0) and Tc{po) hues. For nickel 2so(pNi,0) = 
0.66/0^ and Tcipm) = 627K, so that 7 « 1.25 and 
K « 1.2 X lO^Ka^. 

A rather accurate quantitative look at the itinerant fer- 
romagnetism can be taken via the Stoner approximation 
[3| . In a homogeneous system case, the zero-temperature 
energy per volume is the sum of the kinetic and exchange 
parts: 
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where Pq^ — po/2 ± sq are the up and down spin den- 
sities, J is the Heisenberg exchange interaction con- 
stant [1], and the kinetic energy is defined as IC{pl^) — 

eV(e')de', with e^'^ = ep{pl^) being the kinetic spin 
up and down Fermi energies, and ly being the DOS (per 
spin), [l^l The function ep defines the one-to-one cor- 
respondence between the (spin) density and the position 
of the Fermi energy on the kinetic energy scale. It can 
be defined via its inverse function as p = ^{^')di' , so 
that dtp/dp = v{ep)~'^ and dKL/dp — ep, i.e., ep and p 
are the Legandre conjugates with respect to IC. 

For the ground state, vanishing variations of Sjcx in po 
and So yield the equations determining the ground state's 
spin density and the chemical potential as functions of pq: 



2so(po)^ ^el- e\,, p{po) = [e^p + e\,)/2. 



(3) 



To get (fTa|) from Eqs.dS]) (the graphic solution for pure 
Ni is given in the inset of Figllji) the DOS has to be of 
the following simple form (esd is the 3d band edge): 

iy{e) = i^Ni^'(e3d - e) + i^CuO{e - esd). 

Then, 7 = (z^cu - ^Ni')/(^Ni' + ^^Cu - 2J), and 

p{po) = Mp(Po) - A/i(po), (4) 



where 

Mpo) = esd + lipo- P3d) X I ' ^° < P'"^ (5) 

2 [ ^Cu^ P0> P3d, 

is the chemical potential, corresponding to the paramag- 
netic solution (so — 0, pp ~ p — e^p — ep{po/2)), and 



Am(Po) = (2Az^) ^(pq - P3d)0{p3d ~ Po), 



(6) 



is the chemical potential shift due to the switching 
from the paramagnetic to the ferromagnetic state, with 
Ai^^^ = 7(^7 — i^^i). 14 1 Numerical estimates for 
the parameters of the model can be obtained by fitting 
Mp(PNi,Cu) — Cad from Eq.® and the exchange spitting 
in Ni, 2so(pNi)j7, with their ab-initio values Q: i^q^j^^ = 
6.7,0.79eVa3, J = 1.3eVa3, and Aiy-^ = O.GleVa^.' 

For non-zero temperatures, the free energy of a homo- 
geneous system can be assumed a function of poi sq, and 
T. Its spatial density can be given as the sum of the 
paramagnetic and ferromagnetic (Landau) parts: 

^icx{po,so,T) = Tp{po,T) + AT{po,so,T). (7a) 

J-p is determined mostly by the crystal structure. For not 
very high temperatures, such that all the structural phase 
transitions are far on the (T, p)-plane, J-p can be as- 
sumed temperature independent and consequently equal 
to the paramagnetic part of the zero-temperature energy, 
so that dpgTp = Mp(po)- 

The ferromagnetic properties of the alloy (Eqs.ljT])) de- 
termine the Landau part of the free energy up to an un- 
known factor /: 



A:F = f [2ipo - p3d)'(T - l)(2so/7)' + (2so/7)"j (7b) 

The para-to-ferro shift of the chemical potential provided 
by AF is: 

Apipo, T) = (/(po - P3d)\f - 1)2) 9(1 - f). 

The comparison of the previous equation at T = with 
Eq.®, together with the assumption that / is also tem- 
perature independent, uniquely defines the function /: 



/(po) = (po - P3d)-V(4Ai')- 



(7c) 



Accordingly, Ap(po,T) = Ap(po)(l - T)0(1 - T), with 
Ap{po) from Eq.® (see FigOi). 

Now, in the spirit of the orbital-free local spin den- 
sity approximation (LSDA) we turn to an inhomogeneous 
case by letting the electron and the spin densities spa- 
tially vary: po p{x), so ^ s{x). Capacitively charging 
the conductor we make it an open system, which is gov- 
erned by the fi-functional: 

( Sx{W + Tkx - /ioP + e(p(p - po)/2) . (8a) 
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FIG. 2: (a) The chemical potential provided by the proposed 
model vs. temperature (b) The area on the (T, p)-plane under 
consideration (dashed rectangle). The parameters Ap, pc,5p, 
and 8T are introduced in the text. The arrow symbolizes the 
effect of the (positive) bias on the surface electron subsystem, 
(c) The widths of the spatial profiles of the magnetization and 
the density are determined by the spin and density correlation 
lengths (Eqs[9]), and \a S> Ap. (d) The surface magnetization, 
m, vs. surface charge, a. (insets) The threshold charging, 
CTtr, and the charachteristic magnetization, m*, vs. the alloy's 
detuning from the transition, ST. 



Here Ticx is defined by Eqs.((7]), /io = fi{po,To), 
(p ^ e J (Px'{p{x') — po)l\x — x'\ is the direct inter- 
action potential with e being the electron charge, and 
the nonlocality of the kinetic energy functional is ac- 
counted for by the a- von Weiszacker term Q, W = 
ah"^ / {Sm*)^^^{V p"^ / p]^ , with m* being the effective 

mass of the 3d holes and with p]^ = psd/S — p^^ being 
the spin up and down hole densities. We adopt a = 1/9 
- the case when the a-von Weiszacker approximation is 
asymptotically correct for long wavelengths i.e., the 
domain of the applicability of the Ginzburg-Landau (GL) 
theory below. 

As we already mentioned, the TF layer's spatial elec- 
tron density variation, due to the injection of the carriers, 
is relatively small. Therefore, to be able to reach the fer- 
romagnetic region from an initially paramagnetic state, 
we need an alloy which is close to its transition at Tq: 
Tc = Tcipo) < To, ST^To- Tc « To (see Figllb), so 
that Sp = pq — Pc = ST/ K -c Ap, where pc = Pad — Ap 
and Ap = Tq/k « 0.25/a'^. Having chosen the composi- 
tion, we can focus on a small area on (p, T) plane, such 
that |p — Pel ^ Ap. Within the area, p]j « Ap/2, and 
the von Weiszacker term simplifies as 



ift2((Vp)2 + 4(Vs)2)/(8m*Ap). 



(8b) 



ical potential is given as vtf = fcp'm*/(27r^fi,^), where 
kp — (37r^Ap)^/'^ is the Fermi wave-vector. To make the 
von Weiszacker term energy- wise consistent with the rest 
of the f2-potential, obtained by the physical arguments 
different from the orbital-free LSDA, we can require that 
j^TF = J^Ni, or TO* « 8.2mo. This value of the effective 
mass lies well within the wide range of the 3d sub-bands' 
effective masses, which vary from several mg to almost 
30too (see, e.g., Ref . [lot ) . 

Even though the orbital-free LSDA has proven success- 
ful in some atomic-scale non- homogeneous problems, in 
our case any atomic-scale results obtained from Eqs.([8]) 
would have a rather qualitative character. Indeed, the 
very notion of the composition of a solid solution is well 
defined only on scales larger than the lattice constant. 
Furthermore, the microscopic properties of the electron 
subsystem of the TF-layer differ from those of the bulk 
and are interface-material dependent. The f2-potential 
([5]), however, can be used to find the characteristic dis- 
tances of the spatial variations of p and s, which are 
given by the the density and the spin correlation lengths, 
Ap and As. Omitting intermediate derivations, the cor- 
relation lengths obtained from ([8l) in a linear response 
manner are: 



In the Thomas-Fermi picture of a Fermi liquid, which the 
von Weiszacker correction relies on, the DOS at the chem- 



Xp = 2-i/2ATFReyi + y^l-/ci./(7r2ApA|F), (9a) 
A, = Ap(A/<5r)i/2. (9b) 

Here the Thomas-Fermi radius Atf — (47re^(2z/Ni))^^^^ 
and A = aKj'^Aiy/{4:m*X'j,). 

The numerical estimations lead to A « 7x 10^ K. In the 
vicinity of the transition, the inequality 5T <C A is well 
satisfied, so that the spin and the density scales separate, 
As ^ Ap (see Eq.(|9b|). Thus, we arrived at a typical pic- 
ture of the critical phenomena theory. The spin density 
is a "ready-to-condense" soft variable, behavior of which 
is governed by the large-scale low-energy GL functional 
(the s-dependent part of ([8])), which in the (T, p)-area 
under consideration has the following form: 

Tgl= j (fx{A(Vsf + b{p- pc)s' +Cs^/2), (10) 

with b = 2/{-f'^ApAiy), A = ksb, C = 4V(7^Ap). The 
high-energy stiff variable, p, plays a guiding role via 
the p-dependence of the effective chemical potential for 
the magnetization (the overall coefficient in s^-term of 
Eq.dlOl), B ^ b{p - pc)). The feedback action of the 
soft variable, s, on the stiff variable, p, is weak and/or 
unimportant. 

The proposed form of the near-critical s-p coupling can 
be obtained on a more general footing. Indeed, Tay- 
lor expanding B around the transition point, at which 
B{To,pa) = 0, and noticing that in itinerant ferromag- 
nets the s-p coupling can only be local (especially on the 
As-scale), we arrive at B « b{p — po) -I- b'{To — Tc) = 
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b{p-pc), with PC = po-Sp,Sp = {To-Tc)/K,K = b/b'. 
For quantitative studies, however, it is crucial to possess 
rehable values of the three material-specific parameters 
for Eq. (fTU|) . The way the parameters for Nii_xCux are 
derived in this Letter can now be summarized as follows: 
A is obtained from the orbital-free LSDA considerations 
fEq. ljSbp ): the mutual relation between b and C - from 
the empirical properties of the alloy and from the Lan- 
dau theory for the ILorder phase transitions (Eq.([7b|); 
the overall energy factor for b and C - from previous 
ab-initio parameters of Ni and Cu and from the Stoner 
theory of itinerant ferromagnetism (Eqs.® and fTc]) '). 

On the As-scale, the microscopic effects are scaled out. 
The magnetization is ignorant to the microscopic details 
of the electron density profile, e.g., the Friedel oscilla- 
tions. The injected carriers' spatial density in the TF 
layer, with a width of order of Ap, can be assumed in- 
finitely narrow. Therefore, a positively biased surface 
can be represented as a ID semi-infinite solution (z < 0, 
see FiglJJ:) of the GL Euler equation with p{x) — > po and 
with the following boundary condition imposed by the 
form of the s-p coupling: As S^lnsj^^Q = c/ctr, where 
cr = — {p{z) — po)dz is the surface density of the ex- 
cess holes and ctr = As(Sp is the threshold charging. The 
solution is s = s*(sinh(coth^"'"(o'/(Ttr) — z/Xs))^^8{(J — 
tJtr), where s* = (26(5p/C)^/^. The magnetic response 
of the surface can be characterized by the surface den- 
sity of Borh magnetons, m — j'^ 2s{z)dz (see Figl2}l): 

m = m*log[cr/crtr + ((cr/cTtr)^ - l)^/^]6'(cr - (Jtr), where 
m* = {2A/Cy/^. The characteristic magnetization, m*, 
is independent of 6T, so that the effective magnetic sus- 
ceptibility to the charging, m*/(Ttr oc ST~^^^, grows in- 
finitely as one approaches the transition point (see Figl2ji 
inset). 

In reality, the surface charge is limited by the break- 
down of the insulating interface material. For Si02, the 
breakdown electric field i?si02 ~ lO^V/cm, which corre- 
sponds to the following charging crsi02 = e-E'Si02/(47re) w 
2.5 X lO^'^cm"^, where e = 4.5 is the Si02 dielec- 
tric constant. According to our estimations, for the 
Nii_xCux/Si02 interface, which is ST ~ 5K away from 
the transition, the threshold charging atr ~ 0.15 x 
csi02J on the edge of the semiconductor break- 

down (cr < 173102) the surface magnetization is m w 
1.3 X lO^VBCm-^. 

The magnetic properties of the surface are in fact dif- 
ferent from the bulk on their own, without charging. This 
difference is modeled in the surface phase transition the- 
ory as a delta-functional jump of the local Curie tem- 
perature at the surface of the semi-infinite system (see, 
e.g., Ref.^ and Refs. therein). Within the proposed ap- 
proach the jump can be taken into account as an intrin- 
sic shift of the charging density a CTint + cr (the cases 
Cint ^ are known respectively as ordinary and extrao- 
dinary transitions). A reliable estimate for the material- 



and interface-dependent (Tint can only be obtained from 
ab-initio studies or experimental data. 

Another issue is the quantum fluctuations, due to 
which the critical exponent of the spin correlation length 
must acquire the renormalization group correction: As — > 
XI w Ap(Tc/(5r)V2+'5 ~ 0.14 for, e.g., the d = 4 - e 
approximation [12]). Accordingly, the exponential tail 
of the magnetization profile will be elongated and the 
magnetic response should acquire an enhancement factor 
oc {Tc / . On the other hand, the intrinsic substitu- 
tional irregularity in a solid solution must shorten As due 
to the Anderson localization mechanism. Near the tran- 
sition, the quantum fluctuations' effect is dominant and 
in reality the magnetic response of the surface is stronger 
than the one obtained here on the mean-field level. 

The authors are grateful to Lu J. Sham, Alan H. Mac- 
Donald, and Daniel Neuhauser for useful comments and 
discussions. The work was supported in part by Western 
Institute of Nanoelectronics at UCLA. 



* Electronic address: liovchinn ikov@ ucla.edul 
t Electronic address: |wang@ ee.ucla.edu| 

[1] S. A. Wolf et.al. Science 294, 1488 (2001). 

[2] Y. D. Park et.al. Science 295, 651 (2002); I. Zutic et.al. 
Rev. Mod. Phys. 76, 323 (2004); A. H. MacDonald, P. 
Schifler, and N. Samarth, Nature Mat. 4, 195 (2005); T. 
Jungwirth et.al, Rev. Mod. Phys. 78, 809 (2006). 

[3] H. Ohno et.al. Nature (London) 408, 944 (2000). 

[4] E. C. Stoner, Proc. Roy. Soc. A165,372 (1938). E. P. 
Wohlfarth, Proc. Roy. Soc. A195, 434 (1949). 

[5] S. A. Ahern, M. J. Martin, and W. Sucksmith, Proc. R. 
Soc. (London) 248, 145 (1958). 

[6] L. Michalak, C. M. Canali, and V. G. Benza, Phys. Rev. 
Lett. 97, 096804 (2006); S. Kleff, J. von Delft, M. Desh- 
mukh, and D.C. Ralph, Phys. Rev. B 64, 220401(R) 
(2001); C. M. Canah and A. H. MacDonald, Phys. Rev. 
Lett. 85, 5623 (2000). 

[7] D. A. Papaconstantopoulos, Handbook of the Band Struc- 
ture of Elemental Solids (Plenum, NY, 1986). In the 
estimations, the position of the 3d band edge, esd, 
is taken to be at the energy of the X5 fc-point, exs: 
A'p(pNi,Cu) — esd ~ — 0.21eV and 1.57eV respectively; the 
exchange splitting in Ni is: e^s ~ ^xs ~ 0.82eV. 

[8] R. G. Parr and W. Yang, Density-Functional Theory 
of Atoms and Molecules (Oxford University Press, New 
York, 1989). 

[9] I. V. Ovchinnikov, L. A. Bartell, and D. Neuhauser, J. 
Chem. Phys. 126, 134101 (2007). 

[10] S. Prakash and S. K. Joshi, Phys. Rev. B 2, 915 (1970). 

[11] H. W. Diehl, in Phase Transitions and Critical Phenom- 
ena, vol 10, Eds. by C. Domb and J. L. Lebowitz (Aca- 
demic, London) (1986). 

[12] S.K. Ma, Rev. Mod. Phys. 45, 589 (1973). 

[13] To avoid confusion, even though e]^ > Cp for ferromag- 
netic states, the total kinetic- exchange Fermi energies 
are, of course, the same and equal the chemical potential. 

[14] In fact, the paramagnetic solution always exists. In the 



5 

ferromagnetic region, however, the paramagnetic solution criterion): d^Eicx / dsQ^^ _q < 0. 

is unstable with respect to fluctuations in so (the Stoner 



